Analysis of effects of macroscopic propagation and multiple molecular orbitals on the 
minimum in high-order harmonic generation of aligned CO2 
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We report theoretical calculations on the effect of the multiple orbital contribution in high-order 
harmonic generation (HHG) from aligned CO2 with inclusion of macroscopic propagation of har- 
monic fields in the medium. Our results show very good agreements with recent experiments for 
the dynamics of the minimum in HHG spectra as laser intensity or alignment angle changes. Cal- 
culations are carried out to check how the position of the minimum in HHG spectra depends on the 
degrees of molecular alignment, laser focusing conditions, and the effects of alignment-dependent 
ionization rates of the different molecular orbitals. These analyses help to explain why the minima 
observed in different experiments may vary. 
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I. INTRODUCTION 

High-order harmonic generation (HHG) is one of the 
most interesting nonlinear phenomena when atoms or 
molecules are exposed to an intense laser field [iHl] . Ex- 
perimentally the harmonic fields generated by all atoms 
or molecules within the laser focus co-propagate with 
the fundamental laser field till they reach the detector. 
To compare with experimental measurements, theoreti- 
cal harmonics from each atom or molecule first have to 
be calculated. The induced dipoles are then fed into the 
propagation equation, taking into account the focusing 
condition, the nature of the gas medium, and the equa- 
tion is finally integrated till the harmonics are collected 



HHG from randomly distributed molecules have been 
studied since the 1990's [ol-fl^. In recent years, exper- 
imental HHG studies tend to focus on partially aligned 
molecules (Tsl - flGj . Among the molecules, CO2 is likely 
the most extensively studied system so far [l7l - [20| . Ini- 
tially the interest was focused on the observation of the 
minimum in the HHG spectra of CO2 (2T|-f23|. The posi- 
tions of the minima from different experiments, however, 
are often vastly different. According to the three-step 
model 0, HI] , HHG is generated through the recombi- 
nation of the returning electrons with the molecular ion. 
The interference of electron waves from the two atomic 
centers, under some conditions, results in a minimum in 
the transition dipole, the simplest is the two-center in- 
terference model (26l - l28| . In fact, recombination is an 
inverse process of photoionization. Thus the transition 
dipole is the same as that in photoionization. Any such 
minima have not been observed in photoionization, how- 
ever, since in these measurements molecules are isotrop- 
ically distributed. 

From the theoretical side, the alignment dependence 
of HHG was first studied using the stron g-fi eld approxi- 
mation (SFA) (or the Lewenstein model) [29l - l31j . Subse- 
quently we developed a quantitative rescattering (QRS) 



theory [32M34| for HHG which was a significant improve- 
ment on the Lewenstein model. In QRS the accuracy of 
recombination is treated at the same level as in the pho- 
toionization process. Using QRS, the HHG minimum is 
attributed to the minimum in the photoionization tran- 
sition dipole between the highest occupied molecular or- 
bital (HOMO) and continuum molecular wave functions 
[33, nil . For fixed- in-space CO2 molecules the photoion- 
ization cross sections of HOMO indeed exhibit minima at 
small alignment angles. The experimental HHG spectra 
from partially aligned CO2 have been explained reason- 
ably well by QRS, including the harmonic intensities and 
the phase, so are the polarization and ellipticity of the 
harmonics [H, . 

If HHG is generated from HOMO only, then one ex- 
pects that the position of the minimum does not sig- 
nificantly change with laser intensity according to QRS 
theory. Indeed, strong field ionization depends exponen- 
tially on the ionization potential Ip. The HOMO-1 and 



HOMO-2 orbitals in CO2 
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bound than the HOMO 

to contribute significantly to the HHG spectra 



are 4 and 4.4 eV more deeply 
thus they are not expected 
How- 



ever, it is also well known that tunneling ionization rates 
depend sensitively on the symmetry of the molecular or- 
bital [3^. The HOMO is a ng orbital. It means that at 
small alignment angles the ionization rates are small. For 
HOMO-2, on the other hand, it is a ag orbital, thus it 
has large ionization rate when CO2 molecules are parallel 
aligned. Thus for small alignment angles, HOMO-2 may 
become important even though it is bound 4.4 eV deeper 
than the HOMO. (HOMO-1 is a 7r„ orbital and thus not 
expected to contribute significantly to the HHG.) 

The first step of HHG process is tunneling ionization. 
The alignment dependence of tunneling ionization rate 
is usually calculated using molecular Ammosov-Delone- 
Krainov (MO-ADK) theory [H E^] or SFA. For most 
molecules that have been studied the two models give 
nearly identical alignment dependence (after normaliza- 
tion). However, this is not the case for CO2. Exper- 
imentally, the alignment dependence of CO2 ionization 



2 



reported by Pavicic et al. 4l| is very narrowly peaked 
near alignment angle of 46°. It differs significantly from 
the predictions of MO-ADK and SFA [l^. In fact, so 
far all theoretical attempts [43'-'47'] have not been able 
to confirm the sharp alignment dependence reported in 
the experiment. Furthermore, the observed HHG spectra 
from aligned molecules are inconsistent with the reported 
experimental alignment dependence of ionization 36, 48] . 

The HHG spectra of parallel aligned CO2 have been 
studied in a number of experiments, with 800-nm lasers 
[2H as well as lasers of longer wavelengths [2_3, JQj- It 
was observed that the position of the HHG minimum 
moved to higher photon energies as the laser intensity 
was increased. In Smirnova et al. pTj . such effects were 
interpreted in terms of the multi-channel interference 
(between HOMO and HOMO-2) and the intricate hole 
dynamics. Their calculations were based on HHG gen- 
erated by a single-molecule response. They introduced 
filters into the theory such that only short trajectories 
contributed to the signals. In their calculation, a relative 
phase between different channels due to strong-field ion- 
ization step was introduced "by hand" in order to obtain 
the good agreement with the measurement. 

As noted at the beginning, to compare theoretically 
simulated HHG spectra of molecules with experimen- 
tal ones, the effect of macroscopic propagation should 
be considered. Practically, this has rarely been done. 
For atomic targets, propagation effect is usually included 
with single-atom induced dipoles calculated using the 
Lewenstein model. Only in a few rare occasions the 
atomic response is calculated more accurately by solv- 
ing the time-dependent Schrodinger equation (TDSE) 
[sol [sTI . It is also well-known that the Lewenstein model 
does not predict the HHG spectra (the intensity) pre- 
cisely, but the phases of the harmonics are relatively 
accurate. This fact is used in the QRS model, which 
can be understood as simply replacing the transition 
dipole calculated using plane waves in the Lewenstein 
model by one calculated using accurate molecular contin- 
uum wave functions. The improvement of single-molecule 
HHG spectra calculated using QRS has been well docu- 
mented in our previous publications (3^ . [52j . Compu- 
tationally, QRS is nearly as simple as the Lewenstein 
model (actually QRS is even simpler than SFA for poly- 
atomic molecules, see Ref. [I^), thus induced dipole can 
be easily obtained from QRS to feed into the propaga- 
tion equations to account for medium propagation ef- 
fects. This has been done for atomic targets for low- 
laser-intensity and low-gas-pressure conditions 5J]. It 
has been extended recently to the conditions of high in- 
tensity and high pressure for Ar and to molecular targets 
d, [53, including polyatomic molecules fssj . 

In this paper, we report theoretical studies of the prop- 
agation effect on the HHG of CO2, including contribu- 
tions from HOMO and HOMO-2. Note that experimen- 
tally the effect of propagation on the HHG spectra has 
rarely been investigated either, in particular, its effect 
on the minimum of the HHG spectra. However, this has 



changed recently with the reports of such studies on Ar 
|55l - l57{ . The rest of this paper is arranged as follows. 
In Sec. II, we first summarize the equations used for 
the propagation calculations. We limit ourselves to low 
laser intensity and low pressure. To extend the theory 
to high intensity and high pressure, the optical proper- 
ties of CO2 have to be used and they are not available 
over a broad range of frequencies. We then summarize 
how the single-molecule response from partially aligned 
molecules is calculated. In Sec. HI the calculated re- 
sults are presented. Different factors that can influence 
the precise positions of the HHG minima are examined 
and reported. These results show that precise theoreti- 
cal predictions of the positions of HHG minima in a given 
experiment is difficult, but the trend (the direction of the 
shift of the minimum positions) can be predicted (or ex- 
plained). A short summary at the end concludes this 
paper. 



II. THEORETICAL METHOD 

We first briefiy describe the theory of propagation of 
high harmonics in a macroscopic medium, and in free 
space, till they reach the detector. To calculate the in- 
duced dipole for individual molecules, we include the con- 
tributions from multiple molecular orbitals. Our method 
is based on extending the QRS theory (s^ ]. 



A. Propagation of harmonics in the medium 

In general, both the fundamental laser field and the 
harmonic field are modified when they co-propagate 
through a macroscopic medium. In this paper, we limit 
ourselves to experiments carried out under the conditions 
of low laser intensity and low gas pressure, where the ef- 
fects of dispersion, Kerr nonlinearity, and plasma defo- 
cusing on the fundamental laser field can be neglected 
[1, Hj]. The fundamental laser field is assumed to be a 
Gaussian beam in space. Its spatial and temporal depen- 
dence can be expressed in an analytical form (5^. For 
high harmonics, dispersion and absorption effects from 
the medium are not included. The dispersion and ab- 
sorption coefficients depend linearly on gas pressure and 
could be ignored under low pressure Q . The free-electron 
dispersion is also neglected since the plasma frequency is 
much smaller than the frequencies of high harmonics [H . 

The three-dimensional propagation equation of the 
harmonic field is described as H-ly, N, [541 
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El{r,z',u:,a)=F[El{r,z',t',a)] 



(2) 



3 



and 



Pl{r,z',u:,a)=F[Pl,{r,z',t',a)]. 



(3) 



Here F is the Fourier transform operator acting on the 
temporal coordinate. Eq. ([ij is written in a moving co- 
ordinate frame (z' — z and t' = t — z/c) and the term 

d'^E\ldz''^ is neglected. F^(r-, z\ u, a) is the parallel com- 
ponent of the harmonic field with respect to the polar- 
ization direction of the probe (or generating) laser, and 

pjj(r, z' ,uj, a) is the parallel component of the nonlinear 
polarization caused by the generating laser field, where a 
is pump-probe angle, i.e., the angle between the aligning 
laser and the harmonic generating laser polarizations. 
The nonlinear polarization term can be expressed as 

Pll (r, z',LO, a) = F{[no ~ n,{r, z' , t' , a)]D^^''°\r, z' , t' , a)}, 



where no is the density of neutral molecules, a) 
is the parallel component of the induced single-molecule 
dipole over a number of active electrons [see Eq. (fT2|) 
below], and nc{t',a) is the free-electron density, which 
can be calculated as following: 



nc(t',Q;)= / nc(t',e)p{e, a) sin ede. 



nc{t', 9) = no{l - exp[- 



7(T,0)dT]}, 



(6) 



before detected by the spectrometer. In the meanwhile, 
they may pass through a slit, an iris, or be reflected 
by a mirror. The far-field harmonics can be obtained 
from near-field harmonics through a Hankel transforma- 
tion [eohea] 



El^{rf,Zf,uj,a) = -ik 



E^{r,z',uj,a) krrf 



2{z{ ~ z') 



(7) 



where Jg is the zero-order Bessel function, Zf and rf are 
the coordinates of the far-field points. The wave vector 
k is given by fc = w/c. 

We assume that the harmonics in the far field are col- 
^^■jlected from an extended area. By integrating harmonic 
yields over the area, the power spectrum of the macro- 
scopic harmonics is obtained by 



(5) 



Here Tic{t',6) is the alignment-dependent free-electron 
density, obtained from 



where 7(r, 9) is the alignment-dependent ionization rate, 
which can be calculated by MO-ADK theory [s^, [4^ for 
different molecular orbitals. In Eq. ([5|), 6* is the alignment 
angle with respect to the polarization direction of the 
probe laser, and a) is the alignment distribution with 
pump-probe angle a [H, [sl, HH . 

After the propagation in the medium, we obtain the 
parallel component of near-field harmonics on the exit 
face of the gas jet (z' ~ Zout)- For isotropically dis- 
tributed molecules and partially aligned molecules with 
a = 0° or 90°, by symmetry there are only parallel 
harmonic components. The perpendicular components, 
which are usually much smaller, would appear for par- 
tially aligned molecules and the harmonics will be ellip- 
tically polarized in general [s^ . Generalization of Eq. ^ 
to the perpendicular component is straightforward, but 
we restrict ourselves to parallel component only. Eq. ([1]) 
is solved numerically using a Crank-Nicholson routine for 
each value of w. Typical parameters used in the calcula- 
tions are 300 grid points along the radial direction and 
400 grid points along the longitudinal direction. 



B. Harmonics in the far field 

Once the near-field harmonics are obtained on the exit 
face of the medium, they further propagate in free space 



Sh{uj,a)(x / \El{xf,yf,Zf,uj,a)\'^dx{dyf 



(8) 



where X{ and j/f are the Cartesian coordinates on the 
plane perpendicular to the propagation direction, and 
Tf — \/ x'f + yf. To simulate experimental HHG spectra 
quantitatively, besides laser parameters, detailed infor- 
mation on the experimental setup is needed. 



C. Quantitative rescattering theory for a 
multi-orbital molecular system 

In Eq. laser induced single-molecule dipole mo- 
ment D(t') is calculated quantum mechanically using 
the QRS theory, which has been discussed in detail for 
molecules in Ref. jsj]. Within QRS, laser induced dipole 
moment D{uj,d) for a molecule at a fixed angle 9 (mea- 
sured with respect to laser polarization) is given by 



L>ll-L(tj, 9) = N{9)^/^W{uj)S'^{uj, 9), 



(9) 



where N{9) is the alignment-dependent ionization prob- 
ability, W{uj) is the recombining electron wave packet, 
and d"'^(aj, 9) is the parallel component (or perpendicu- 
lar component) of the photorecombination (PR) transi- 
tion dipole (complex in general). Here we only consider 
linearly polarized lasers and linear molecules. W{u!) is 
independent of the alignment angle 9. From Eq. (9), it 
can be expressed as 



W{uj) = 



£)II-L(a;,6l) 
N{9y/^dl^{uj,0)' 



(10) 



In QRS, W{uj) is usually calculated only once for a 
given angle 9 using SEA, where all the elements on the 
right-hand side of Eq. ([T0| are obtained by replacing the 
continuum waves with plane waves. In QRS, accurate 
c?"'^(aj,0) are obtained from quantum chemistry code 
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[63, [6J] and tunneling ionization probability N{9) are ob- 
tained either from MO-ADK or from SFA. Put all 
of these together into Eq. ([9]), laser induced dipole mo- 
ment D{uj, 6) for each orbital is obtained. Note that in 
QRS, the wave packet W{ijj) automatically includes the 
phase which is dependent on the molecular orbital. Thus 
there is no need to introduce the relative phase between 
orbitals, in contrast to the approach in Ref. [2l|. SFA is 
used to calculate the ionization probability in Eq. ([9]) in 
this paper unless otherwise stated. 

We assume that the degree of molecular alignment 
is not varied spatially within the medium because the 
molecules are usually partially aligned by a weak and 
loosely focused laser beam [HI] . By coherently averaging 
the induced dipole moments over the molecular angu- 
lar distribution, we obtain the averaged induced dipole 
of partially aligned molecules at each point inside the 
medium: 

/■TT 

D^^'^^^iuj^a)^ D\\u:,e)p{e,a)smede, (11) 

where p{9, a) is again the angular (or alignment) distri- 
bution of the molecules with respect to the polarization 
direction of the probe laser. For randomly distributed 
molecules, p{9,a) is a constant. Note that the above 
procedure is only for the specified molecular orbital. 

For the specific problem addressed in this paper, we 
consider electrons either in the HOMO (Itt^) or in the 
HOMO-2 (3cr„) of CO2. The electrons are freed and 
then recombine back to the same orbital. As discussed 
above, the multiple orbital effects are important at small 
alignment angles only due to symmetry consideration. 
At these angles HOMO-1 (ItTu) is ignored since the ion- 
ization rate of HOMO-1 is quite small compared with 
HOMO and HOMO-2 [H Eo, |4|. At large alignment 
angles, only HOMO becomes important. The total laser 
induced dipole over a number of active electrons can be 
written as [69, rool 



are fed into Eg. ^ 
Jin et al. 



The same procedure is used in 
54| for atomic targets. Note that in the 



7?ll^*°*(w,a) = 



(12) 



where index j refers to the HOMO and HOMO-2 of CO2, 
and n is an index to account for degeneracy in each molec- 
ular orbital. 

Before proceeding, we comment that the QRS theory is 
formulated in the energy (or frequency) domain. There is 
no explicit treatment of "core dynamics" as in Ref. [2l| . 
The time evolution of the core is reflected only by the 
energy of the core in the free propagation of the electron 
wave packet. In other words, any possible many-body 
interchannel couplings between HOMO and HOMO-2 are 
not included in the present treatment. Before such effects 
are addressed, other factors that are more important on 
the HHG as discussed here have to be treated first. 

For the propagation of harmonics in the medium, we 
need to obtain hundreds of the total induced dipoles 
L)ll'*°*(a;, a) for different laser intensities, and then they 



present model, dielectric properties of molecules due to 
non-isotropic distributions are also neglected. 



III. RESULTS AND DISCUSSION 

A. HHG spectra of randomly distributed CO2: 
theory vs experiment 

HHG spectra by 800-nm and 1200-nm lasers have 
been reported for isotropically distributed and partially 
aligned N2 and CO2 molecules (43. The spectra of N2 
have been analyzed by Jin et al. i8!,l55| recently including 
only the HOMO. 

In Figs. [Ha)-[TJc), we show the HHG spectra for 
isotropically distributed CO2 molecules by 800-nm and 
1200-nm lasers. To obtain good agreement between the- 
ory and experiment, especially in the cutoff region, the 
intensity used in the theory is adjusted from the value 
given in the experiment. In Figs. [Ha)-[TJc), the inten- 
sities in theory (experiment) are 1.9 (2.1), 0.8 (1.0), 1.0 
(1.2), in units of lO^^W/cm^, respectively. Other param- 
eters used in the simulation are the same as those given 
in the experiment [i^. The laser parameters are: pulse 
duration is ^ 32 fs (800 nm) or ^ 44 fs (1200 nm), beam 
waist at the focus is ^ 40 pm. A 0.6-mm-wide gas jet is 
located 3 mm (800 nm) or 3.5 mm (1200 nm) after the 
laser focus, and a slit with a width of 100 pm is placed 
at 24 cm after the gas jet. 

Figs.[TIa)-[IJc) clearly show the good overall agreement 
between experiment and theory for randomly distributed 
CO2 molecules. We have checked that HOMO is domi- 
nant for randomly distributed CO2, with negligible con- 
tributions from inner orbitals. The HHG spectra do not 
exhibit any minima, as opposed to the spectra of random 
N2 molecules when they are generated under the same 
experimental conditions [49]. For randomly distributed 
CO2, there was no minimum found in HHG spectra us- 
ing an 800-nm laser by Vozzi et al. :22j. However, for 
the 1300-nm lasers, the data from Torres et al. [2^ H^l 
appear to show a very weak minimum at photon energy 
near 45 eV. Without more careful study including dif- 
ferent intensities and wavelengths, however, this is not 
conclusive. 



B. HHG spectra of aligned CO2: theory vs 
experiment 

Experimentally HHG spectra have also be reported 
from aligned CO2 molecules. A relatively weak and 
short laser pulse was used to impulsively align molecules, 
and the HHG spectra were taken at half-revival (~ 21.2 
ps in CO2) when the molecules were maximally aligned 
[49| . The angular distributions of the aligned molecules 
are obtained by solving the TDSE of rotational wave 
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Photon energy (eV) 

FIG. 1: (Color online) Comparison of theoretical (green lines) and experimental (red lines) HHG spectra of random and 
aligned CO2 molecules, in an 800-nm laser shown in (a) and (d), and in a 1200-nm laser shown in (b), (c), (e), and (f). Laser 
intensities are indicated where Io=10^'* W /ere? . Experimental data are from Ref. [4^. Arrows indicate the positions of minima. 
Pump-probe angle a=0°. See text for additional laser parameters and experimental arrangements. 



packet nil. The degree of alignment is (cos^0)=O.6O in 
Fig.lUd), and (cos2 6')=0.50 in Figs, [^e) andlllf). The 
polarizations of the pump and probe lasers are parallel 
to each other. 

The HHG spectra of partially aligned CO2 molecules 
are shown in Figs. [lld)-[ljf), which are obtained under 
the same generating lasers and experimental arrange- 
ments as those in Figs.[IJa)-[TJc), respectively. The simu- 
lation and experimental data agree well with each other 
in general. In Fig.[lje), the discrepancy is a little bigger, 
showing the drop near 40 eV is larger in the experiment 
than in the theory. But we note that in Fig. [Hf), the 
experimental data do not drop as rapidly, in agreement 
with the theoretical simulation. 

The minima in the HHG spectra of CO2 and their de- 
pendence on laser intensity have been widely discussed 
in the literature [H ill. In Fig. [J^d), for an 800-nm 
laser, experiment gives a strong minimum at 42±2 eV, 
our simulation predicts a minimum around 42 eV. For 
the 1200-nm laser, in Fig. [Ije), experiment shows a min- 
imum at 51 ±2 eV, theory predicts a minimum around 
50 eV. In Fig. [Ijf), the experimental minimum is shifted 
to 57±2 eV, and the theoretical one is moved to around 
53.5 eV. Thus our simulation also shows the shift of the 
minimum from low to high harmonic orders as laser in- 
tensity is increased. Below we interpret the origin of the 
shift. 



C. Origin of minimum in the HHG spectra of 
aligned CO2: Type I and Type II 

We next analyze the origin of the minimum in the HHG 
spectra seen in Figs.[ljd)-[ljf), and consider the dominant 
contributions from HOMO and HOMO-2 only. First we 
define the averaged PR transition dipole for each molec- 
ular orbital by 

dii'^^s(^^Qf) = / N[eY''^S[uj,e)p{e,a)sinede. (13) 

Jo 

This gives a measure of the relative contribution of 
each molecular orbital to the HHG, which is obtained by 
averaging over the angular (or alignment) distribution of 
the partially aligned molecules, weighted by the square 
root of the tunneling ionization probability. The relative 
ionization rates between HOMO and HOMO-2 change 
with laser intensity. 

Figs. [lla)-mc) show the envelopes of the HHG spectra 
from individual molecular orbitals together with the total 
ones, each obtained after propagation in the medium. 
In the meanwhile, the averaged PR transition dipoles of 
HOMO and HOMO-2 under different generating lasers 
and alignment distributions are shown in Figs.[2jd)-[2jf), 
respectively. 

In Figs. [2ja) and [2jb), there are no minima in the 
HHG spectra of HOMO or HOMO-2, but the minimum 
shows up in the total spectra. This is due to the in- 
terference between HOMO and HOMO-2. We call this 
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FIG. 2; (Color online) (a)-(c) Macroscopic HHG spectra (envelope only) corresponding to Figs. [TJd)-[TJf), respectively. Total 
(HOMO and HOMO-2 together) spectra and the spectra of individual HOMO and HOMO-2 are shown, (d)-(f) Averaged 
photorecombination transition dipoles (parallel component, the square of magnitude) of HOMO and HOMO-2 corresponding 
to (a), (b), and (c), respectively. Laser intensities are indicated where Io = 10^* W/cm^. Arrows indicate the positions of 
minima. Pump-probe angle a=0°. 



type I minimum. Clearly the minimum position will 
change with laser intensity since the relative ionization 
rates between HOMO and HOMO-2 change with inten- 
sity [also see Figs.jSfc) and[SUd)]. Similar analysis can be 
found in Ref. 2l[. In Fig. [2][c), there is a minimum in 
the HOMO spectra at 52.6 eV. This minimum is shifted 
to 53.6 eV in the total spectra due to the interference 
with the HOMO-2. This is categorized as type II mini- 
mum. Similar analysis of this type can be found in Refs. 
[23I . [68j . The minimum in the HOMO spectra is due 
to the minimum in the averaged PR transition dipole of 
HOMO shown in Fig.[2If). But their positions differ due 
to modification of the macroscopic wave packet (MWP). 
In this connection we mention that the earlier calcula- 
tions [3^ . Issj with an 800-nm laser showed minimum in 
HHG spectra at small pump-probe angles due to the con- 
tribution from the HOMO only. These calculations were 
carried out with a higher degree of alignment and higher 
laser intensities as compared to our present study. This 
is expected as the minimum in the averaged PR tran- 
sition dipole from HOMO in Fig. ^d) becomes deeper 
and is slightly shifted away from the cutoff to lower ener- 
gies with increased degree of alignment [see Fig. [5Je) and 
mf)]. Furthermore, as shown in Ref. |55| . the minimum 
in the HOMO spectra could also result from the multipli- 
cation of MWP and averaged PR transition dipole even 
when neither has minimum. When a minimum occurs in 
the dominant orbital, its position will not change much 
if it remains the dominant one when the laser intensity 



changes. The little bump around 36 eV in the HOMO 
spectra as well as in the total spectra can be seen due to 
the bump in the HOMO curves in Figs.[i;d)-[i;f). Its po- 
sition does not change much since the HOMO-2 remains 
small. 



Previously in Refs. [5a], we have shown that the 
macroscopic HHG spectrum is the product of a MWP 
and an averaged PR transition dipole for each individual 
molecular orbital. Since the ionization rate for each or- 
bital has been incorporated in the averaged PR transition 
dipole, the MWP is mostly identical except for the phase 
due to ionization potential. The averaged PR transition 
dipole is very sensitive to ionization rates. The relative 
magnitude changes rapidly with the increase of laser in- 
tensity. Thus when two averaged PR transition dipoles 
are comparable [see Fig. Hfd)], the position of the mini- 
mum changes rapidly with the laser intensity. The aver- 
aged PR transition dipole is also sensitive to alignment 
distributions [see Figs. [2ld)-[2lf)]. At low laser intensity, 
HOMO-2 is small, interference often occurs in a narrow 
region only where the two amplitudes are comparable, 
see Figs. IDJb) and[2{c). In comparison, in Smirnova et 
al. HOMO and HOMO-2 tend to interfere over a 

broad photon-energy region. The ionization rates and 
transition dipoles used in their calculations are different 
from ours. 
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FIG. 3: (Color online) Laser intensity dependence of macro- 
scopic HHG spectra (envelope only) (a) in an 800-nm laser 
and (b) in a 1200-nm laser. Intensities are shown in units 
of lo^lO^* W/cm^. Arrows indicate the positions of minima. 
Degree of alignment: (cos^ 6)=0.60. Pump-probe angle q:=0° . 



D. Progression of harmonic minimum vs laser 
intensity 

In Figs. EJa) and Elb) we show the envelope of the 
calculated HHG spectra for four different peak intensi- 
ties with an 800-nni laser and a 1200-nni laser, respec- 
tively. For the 800-nin spectra, the lowest intensity does 
not have a minimum. For the higher ones, each spec- 
trum has a type I minimum, with its position shifts to 
lower photon energy as the laser intensity is decreased. 
The degree of alignment of molecules used in the cal- 
culation is (cos^ 6')=0.60. We find that the shift can- 
not be attributed to either MWP or the averaged PR 
transition dipole alone. For the 1200-nm data, also with 
(cos^ 9)— 0.60, which is different from Figs.[TJe) and[IJf), 
we find that the minimum is type II, where the averaged 
PR transition dipole of the HOMO has a minimum. The 
minimum in the HHG spectra of the HOMO shifts to 
higher photon energy as the intensity increases, but the 
interference with HOMO-2 shifts the minimum to even 
higher energies. In other words, the shift of the position 
of the HHG minimum vs intensity cannot be attributed 
to a single factor alone. 



E. Other factors that influence the precise 
positions of HHG minima 

In our analysis, the averaged PR transition dipole is 
calculated over the angular distribution of the molecules 
and thus depend on the degrees of alignment. Since the 
latter cannot be accurately measured, we check how sen- 
sitive the calculated spectra is with respect to the as- 



sumed alignment distribution. In Fig.|4jc), four different 
alignment distributions are shown. The distributions are 
multiplied by the volume element sin 9 for easy compar- 
ison. Three of them are obtained from the calculated 
rotational wave packets [58], with {cos^ 6) as 0.63, 0.60, 
and 0.55, respectively. The other is the commonly used 
cos"* 9 distribution. For 800-nm and 1200-nm lasers, the 
envelopes of the calculated HHG spectra are shown in 
Figs.m^a) andlUJb), respectively. The precise position of 
the minimum changes slightly except for the one from 
the cos^ 9 distributions. However, change of a couple of 
eV's is seen. 

To precisely determine the minimum in the HHG spec- 
tra, accurate alignment- dependent ionization probability 
N{9) for each molecular orbital is needed. For CO2 ■ even 
for HOMO, different theories in the literature [21I 113- 
Sl, I45l - l47| show non- negligible differences, and they do 
not agree with the experimental data . Here we exam- 
ine how the HHG spectra change with the different ion- 
ization rates used. The ionization rates for both HOMO 
and HOMO-2 can all be easily calculated from SFA or 
from MO-ADK theory. Figs. EJa) and Efb) show the 
HHG spectra calculated using the ionization probabili- 
ties shown in Figs.[5l[c) and[5jd). Other laser parameters 
used in the calculation are given in the figure captions. 
The difference of the position of the minimum is 3 eV in 
Fig.EJa) and 2 eV in Fig. [Sl^b). Note that the ionization 
probabilities from SFA and MO-ADK are normalized at 
the peak of the HOMO curve. In Fig. [5ja) the spectra 
are normalized at H33 (51 eV) and in Fig. EJb) at H65 
(67 eV). 

The HHG spectra are also sensitive to the experimental 
arrangement and thus can also move the position of the 
HHG minimum. To demonstrate this, we (i) move the 
gas jet to the laser focus and collect the signal using a 
slit; (ii) put the gas jet after the laser focus, and collect 
HHG signal without slit (total signal). These two will 
be compared to the arrangement used in this paper: gas 
jet is after the laser focus and the HHG is collected with 
a slit. The results are compared in Fig. [6l Note that 
the spectra are normalized at H17 (26 eV) in Fig. [6lja) 
and at H35 (36 eV) in Fig. MJ^)- Not only the spectra 
change quite significantly, but also the position of the 
HHG minimum. This illustrates that it is very difficult to 
compare the position of the HHG minimum from different 
experiments. In this comparison, the change of the HHG 
spectra is due to the change of MWP which depends 
on the experimental setup. The averaged PR transition 
dipoles of HOMO and HOMO-2 are the same in the three 
calculations. 



F. The dependence of the HHG minimum on the 
pump-probe angle 

In Worncr et al. [49] , it was found that the minimum 
in the HHG spectra of aligned CO2 moved to higher pho- 
ton energy with increasing pump-probe angle a, i.e., the 
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FIG. 4: (Color online) Dependence of macroscopic HHG spectra (envelope only) with degrees of molecular alignment distribu- 
tions for (a) an 800-nm laser with intensity of 1.8x10^'' W/cm^, and (b) a 1200-nm laser with intensity of 1.0x10^^ W/cm^. 
The weighted angular distributions of the molecules are shown in (c). Arrows indicate the positions of minima. Pump- probe 
angle q=0°. 
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FIG. 5: (Color online) Dependence of macroscopic HHG spectra (envelope only) on the ionization probabilities calculated from 
MO-ADK or SFA in (a) an 800-nm laser, and (b) a 1200-nm laser. Laser intensities are indicated where 10 = 10^** W/cm^. 
Arrows indicate the positions of minima. Degree of alignment: (cos'^ &)=0.60. Pump-probe angle a=0° . (c) and (d) Alignment- 
dependent ionization probabilities of HOMO and HOMO-2 calculated using MO-ADK and SFA. Laser parameters are the same 
as (a) and (b). Ionization probabilities of HOMO-2 in (d) are multiplied by 5. 



angle between aligning pump beam and the HHG gener- 
ating probe beam polarizations. This phenomenon has 
also been observed in other measurements j2l], |6^. In 
contrast, the minimum in the HHG spectra of aligned 
N2 has been found to stay at the same position for all 
pump-probe angles [H, [69| . 

We show the calculated HHG spectra at four pump- 
probe angles in an 800-nm laser in Fig. [71[a), and in a 
1200-nm laser in Figs. [71[b) and[71[c). Laser parameters 
and experimental arrangements are the same as in Fig. [T] 
and the degree of alignment is (cos^ 6')=0.60. The HHG 
spectra for a=0° have been shown in Fig. O These fig- 



ures show that the minimum in the HHG spectra moves 
to higher photon energies with increasing a , and the 
minimum disappears at large angles. At larger pump- 
probe angles, contributions from large alignment angles 
increase. Since HOMO dominates over HOMO-2 at large 
angles in both the ionization rates [see Figs. [Sfc) and 
[SJd)] and the PR transition dipoles [see Figs. [7Ud)-[7Uf)], 
thus HHG at large pump-probe angles has essentially no 
contributions from HOMO-2. Also note that at small a, 
the total HHG yield is much smaller [sj]- In fact, the to- 
tal HHG spectra for randomly distributed CO2 have little 
contributions from molecules that are aligned nearly par- 
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FIG. 6: (Color online) Dependence of macroscopic HHG spec- 
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dashed lines). Arrows indicate the positions of minima. De- 
gree of alignment: (cos^ 6')=0.60. Pump-probe angle a=0°. 



allel to the polarization axis of the laser. The HHG spec- 
tra of CO2 are complex only in the region where HHG 
yields are small. This is generally true - interpretation 
of small processes always requires careful and detailed 
theories. 



IV. SUMMARY AND OUTLOOK 

In this paper we have analyzed the multiple orbital 
contribution to HHG in CO2 with the inclusion of macro- 
scopic propagation effect. In the past few years, there 
have been many experimental and theoretical studies 
on the HHG of CO2 from many laboratories, using 
lasers with different wavelengths and intensities, for CO2 
molecules that are randomly distributed or partially 
aligned. In particular, for CO2 molecules that are par- 
tially aligned along the polarization axis of the probe 
laser, many experiments have shown that the HHG spec- 
tra exhibit minima and the positions of the minima shift 
with laser intensities [2l|, [2I, . The shift of the mini- 
mum position with laser intensities has been attributed 
to the interference between the contributions to HHG 
from the HOMO-2 with the one from the HOMO, despite 



HOMO-2 lying at 4.4 eV deeper than the HOMO. Since 
HHG is a nonlinear process, these observations posed a 
great challenge to the theory, especially for the predic- 
tion of the position of the minimum and how it changes 
with the experimental conditions. Since all experimen- 
tal HHG spectra include macroscopic propagation effect, 
comparison of theory with experiment is incomplete un- 
less the theory also has included the propagation effect. 
Our analysis in this paper is based on the macroscopic 
propagation code extended for aligned molecules, and the 
recently developed quantitative rescattering theory. We 
find that although HHG spectra change significantly un- 
der different experimental parameters such as degree of 
alignment, focussing condition, and the use of a slit, the 
position of the minimum in the HHG spectra behaves in 
a similar trend as laser intensity and pump-probe angle 
vary. This trend has been found to be consistent with the 
recent experimental measurements from different groups. 

We comment that the present theory and the earlier 
one by Smirnova et al. \2l\ both explain the intensity 
dependence of the change of HHG minima, but the de- 
tails between the two theories are quite different. The 
alignment dependence of the ionization rates, the recom- 
bination dipole matrix elements and their phases entering 
the two theories are not the same, for both the HOMO 
and HOMO-2. As illustrated in this paper, these param- 
eters can all affect the position of the predicted inter- 
ference minimum. Furthermore, in Smirnova et al. (2]| 
the interference is attributed to the importance of hole 
dynamics in the ion core. Our approach is formulated 
in the time-independent fashion, no hole dynamics is in- 
cluded. Since HHG spectra are taken without explicit 
observation of electron dynamics, the difference between 
the two models cannot be settled. Despite these differ- 
ences, our understanding of the HHG spectra of CO2 has 
come a long way since 2005 [1^. With the possibility 
of including macroscopic propagation effect "routinely" 
in the HHG theory for molecular targets, further experi- 
mental studies should explore the effects of laser focusing 
condition and gas pressure, for lasers extending to longer 
wavelengths. Such studies would further our basic under- 
standing of strong-field physics of molecules to the next 
level, and eliminate the need of introducing extraneous 
assumptions in the interpretation of experimental data. 
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